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Abstract 

We study the formation and evolution of a turbulent spectrum of Alfven waves driven by reflection off the 
solar wind density gradients, starting from the coronal base up to 17 solar radii, well beyond the Alfvenic 
critical point. The background solar wind is assigned and 2D shell models are used to describe nonlinear 
interactions. We find that the turbulent spectra are influenced by the nature of reflected waves. Close to the 
base, these give rise to a flatter and steeper spectrum for the outgoing and reflected waves respectively. At 
higher heliocentric distance both spectra evolve toward an asymptotic Kolmogorov spectrum. The turbulent 
dissipation is found to account for at least half of the heating required to sustain the background imposed solar 
wind and its shape is found to be determined by the reflection-determined turbulent heating below 1 .5 solar 
radii. Therefore reflection and reflection-driven turbulence are shown to play a key role in the acceleration of 
the fast solar wind and origin of the turbulent spectrum found at 0.3 AU in the heliosphere. 

Subject headings: MHD — waves — turbulence — solar wind 



1. INTRODUCTION 

Recent high resolu tion observations from Hinode 
(iDe Pontieu ei al.l l2007h and the Coronal Multi-channel 
Polarimeter (iTomczyk et al.l l2007h have shown that the 
solar atmosphere is pervaded by Alfvenic (or kink-like, e.g. 
|Van Doorsselaere et al. 2008) oscillations: observed in jets, 
spicules or in coronal loops, velocity and magnetic field 
oscillations {Sb, 6u) are coupled and propagate at speeds 
close to the Alfven speed. Far from the sun, between 0.3 AU 
and several AU, in situ data show that in the frequency range 
10 Hz < / < 10"^ Hz, fluctuations in magnetc field and 
velocity 6b and 6u possess many of the properties of outward- 
propagating "spherically polarized" Alfven waves, namely: 
quasi incompressibility, correlated oscillations, and a con- 
stant (total) magnetic field intensity, while at the same time 
revealing their turbulent nature through a well-developed 
power-law frequency spectrum, with a break separating 
different power law slopes of -1 and -1.6 which moves 
to lower frequencies with increasin g distance from the sun 
(iBavassano et al.ll982l:lfu et al.ll98 4). The 5m -5^7 correlation, 
upon which the propagation direction determination is made, 
depends on the frequency considered an d varies w ith distance 
( IBavassan o et al. 2000b a) and latitude (iGrappinl l2002). typi- 
caUy in flie range 1 /2 < \(6u ■ 6b I ^jATip)l{6u^ + 6b^lAnp)\ < 1 
(perfect correlation corresponding to 1). These facts suggest 
that the inward travelling wave-mode component, required 
for nonlinear couplings between incompressible fluctuations, 
must indeed be present. 

This component might be generated locally between 0.3 AU 
and 1 AU by shear, compressible or pick-up ions interactions, 
or it could be already formed in the sub-Alfvenic corona and 
later on nonlinearly advected into the heliosphere by the solar 
wind, the hypothesis we consider here. 

' Observatoire Royale de Belgique, 3 Avenue Circulaire. 1180, Bruxelles, 
Belgium; e-mail: verdini@oma.be 

- Dipartimento di Astronomia e Scienza dello Spazio, Univ. di Firenze, 
Largo E. Fermi 3, 50125, Firenze, Italy 

' Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak 
Grove Drive. Pasadena, CA 91 109. USA 

Institut d'Astrophysique Spatiale, CNRS - Universite Paris Sud, Bat. 
121, 91405, Orsay Cedex, France 



The dynamics inside the Alfvenic point region is of primary 
importance to understand the origin of the spectrum one finds 
at 0.3 AU and whether it has any role in accelerating the 
solar wind. The variation of the propagation speed induced 
by density gradients in the stratified corona and accelerating 
solar wind causes outward Alfven waves to be reflected, 
predominantly at lower frequencies, hence triggering the 
incompressible cascade. The power dissipated by the cascade 
contributes to coronal heating, also modifying the overall 
turbulent pressure gradient, fundamental to the acceleration 
of the fast solar wind. 

While there are several studies on the linear propagation 
and reflection of Alfven wa ves in the sub-A lfvenic coron a and 
solar wind (.Hollweg 19781 l Hememami & Olbert 19801: IWlil 
IT99I iHollweg & Isen berg 20017 excepting phenomenolog- 
ical models with an essentially dimensional estimate of the 



role o f turbulent heating ((H ollweg et al. 1982; Dmitru k et al 
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2007h very few of them have considered nonlinear interac- 
tions. IVelh et al.1 (Il989l Il990h studied the turbulent cascade 
sustai ned by reflected wave s in the super- Alfvenic solar wind, 
while iDmitruk et al.l (|2002|) considered the same mechanism 
in the sub-Alfvenic corona below 3 Rq, hence neglecting the 
solar wind. 

In the present letter we will extend these studies follow- 
ing the development of the turbulent cascade from the base 
of the corona up to 17 Rq, well beyond the Alfvenic crit- 
ical point (located at about 13 Rq in the solar wind model 
adopted). Direct numerical simulations are still prohibitively 
costly in terms of computational times, so n onlinear interac- 
tions are simulated using a 2D shell model ( Buchlin & Vellil 
l2007h which simplifies nonlinear interactions but still allows 
4 decades in the perpendicular wavenumber space to be cov- 
ered while rigorously treating the propagation and reflection 
of waves along the radial mean magnetic field. 

2. MODEL DESCRIPTION 

The equations describing the propagation of Alfven waves 
in an inhomogeneous medium are derived from magnetohy- 
drodynamics (MHD), assuming that the large scale fields are 
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stationary and separ ating the time-fluctuating fields from th e 
large-scale averages (iHeinemann & Olbertlll980l:lwiil[T993h . 
Therefore the large scale magnetic field, bulk wind flow and 
density (B, U, p respectively) appear as specified coeflicients 
in the MHD equations for the fluctuations. 

We consider a magnetic flux tube centered in a po- 
lar coronal hole, which expands super-radially with a 

i non-dimensional) area A(r) - r^f(r) first p rescribed by 
Kopp & HolzHlll976t iMunro & JacksonlllQTTl) . Distances 
are normalized to the solar radius, and the coefficients are 
given by ro = 1.31, cr - 0.5, f,„a.x = 7.26 (respectively the 
location, width, and asymptotic value of the super-radial ex- 
pansion) so that A(l) =1. The field becomes B(r) = B^/Air) 
where we take Bq = 10 G. The wind speed and density 
U{r), p(r) are obtained solving the ID momentum equation 
with an assigned temperature T{r), Tq - 4 10^ K and a nu- 
merical density at th e coronal base Hq - 5 10^ cm"^ (see 
IVerdini & Vellill2007l and references therein for details on the 
equation and on the temperature profile). The resulting wind 
is supersonic far from the sun with Uiau ~ 750 km/s and 
"lAU ~ 3 cm"^. The Alfven critical point (ra) is at about 
13 Rq, the sonic critical point is at about 1.9 Rq, the Alfven 
speed Va - B/ -\j4np at the base is Va^o ~ 1000 km s"' and 
has a maximum Va,,na\- ~ 3500 km/s at r = 1.5 Rq- At the end 
of the domain (17 Rq) U ~ 740 km/s and Va ~ 630 km/s. 
We assume 6u to be incompressible and transverse with re- 
spect to B. The momentum and induction equations for 
6u and 6b are written in terms of the Elsasser fields Z'^ - 
6u + 6b/ ^jAnp, corresponding to Alfven waves which propa- 
gate respectively outwards and inwards in the solar wind ref- 
erence frame. 

Substituting the nonlinear terms, which act in planes per- 
pendicular to the radial direction, wit h a 2D MHD she ll model 
representation (in the form given by iBiskampll 19941) and as- 
suming radial propagation finally yields the model equations 
forZ^ =Z±(r,;tJ =Z,^(r): 

The complex scalar values, m„ - (Z,| + Z,7)/2, b„ - (Z,7 - 
Z,^)/2, represent the velocity and magnetic (in velocity units) 
field fluctuations corresponding to the scale An - AqI^" - 
2n/k„, n is the shell index, and accounts for nonlinear in- 
teractions of the form Z'^Z'^^, with l,m - n + \,n ±2. Finally, 
= (v + ?7)/2 are combinations of the kinematic viscosity 
and the magnetic resistivity (we take rj - v). 
Simulations are carried out with the code Shell-Atm^ 
jBuchlin & Velli 2007). The advection terms in eq.[T]are com- 
puted with a second order upwind scheme (Fromm scheme) 
which allows a good conservation of the phase of the fluc- 
tuation. Time is advanced with a third order Runge-Kutta 
for the nonlinear part of the equations. The radial domain 
is decomposed in ~ 25,000 planes over a non-uniform grid 
while 21 shells are used for the nonlinear interactions. Trans- 
parent boundary conditions are imposed at the top for both 
waves, and at the bottom for the Z . Here all the gradi- 
ents are artificially set to zero, in order to avoid reflections. 
Energy is injected in the domain imposing the amplitude 
- fnit) the first 3 shells corresponding to length scale 

' modified to include the wind spherical expansion 



of the order of 8.000-34.000 km (/lo.o = 0.02 Rq in the shefl 
model), with /(f) a function with a time correlation and pe- 
riodic ity r* =a 1000 s. The form of the function is given in 
Buchli n & Vellil ( l2007h . here it is important to note that de- 
spite the correlation time is only 1000 s, some low-frequency 
fluctuations are injected for long time series (as in the present 
simulation). Simulations last about 20 crossing time scale, 

Tcr - ~ 7000 s, and time averaged quantities are com- 

puted on the last 10 Tcr, during which the system has an ap- 
proximate stationary state. 

3. REFLECTION AND NONLINEARITIES 

Given that we start with an initial flux of Z^ at the coro- 
nal base, reflection is the only trigger for nonlinear interac- 
tions. In Figure [T] the root mean squared (rms) amplitudes 

l-Z'^l - VSn l-^jfP plotted as a function of time and he- 
liocentric distance. Two components of the reflected waves 
are clearly visible in the Z contours: a "classical" compo- 
nent Z^^^j^, which propagate with the expected U - Va phase 
speed (negative below the Alfvenic critical point, ra » 13 Rq) 
and an "anomalous" component Z^^om^ which travels with the 
same speed as the mot her wave, U + Vg. As shown analyti- 
cally and numerically dVelh et al.lll990l : [Hollweg & IsenbergI 
l2007h this "anomalous" component is the direct product of re- 
flection. Generally in presence of density gradients, for small 
values of the ratio e = a/oj, typical of the upper corona, with 
a ^ Tf} ^ \(U + Va)V'alVa\ in eq. [T] each field can be de- 
composed in a primary and secondary component. A Z"^ pri- 
mary component is injected at the base, while the Z is made 
up of only the secondary component, given by reflection. In 
each plane, as the Z^ arrives, the secondary component can 
be seen as the result of a forcing term given by » aZ^ , hence 
producing a wave which travels with the same phase speed 
and of the above amplitude. The value Z - eZ^ follows nat- 
urally by finding the "forced" solution to the linearized equa- 
tion for Z . At later time, as the Z^ has propagated away, the 
forcing disappears and the secondary component propagates 
backward with the classical phase speed. When a Z^ pulse 
wave is excited in the corona, Z appears as a halo spreading 
backward from the mother Z^ wave. Nonlinear interactions 
modifies the above picture, acting as a local (in a given plane) 
source for the Z,7 which is uncorrelated with respect to that 
given by reflection and hence generating waves propagating 
with classical phase speed U - Va, i.e., a primary component. 
Inward propagating waves have a very long propagation time 
at ra that slows down the overall relaxation toward a steady 
state: for example in Figure [1] the increase of the Z ampli- 
tude at f ^ %Tcr results from the superposition of the Zg^om 
with the backward propagating Z^,^^^ produced at f « 3t„ . 

The different nature of the reflected waves influences the 
spectral energy transfer. Generally speaking, while the wave 
amplitudes increase with distance the nonlinear timescales 
~ (^n-Z,^)"' decrease, following the reflection coefficient 
and the flux tube expansion. At small wavenumbers r^^ < r^,, 
the energy transfer is more efficient for the reflected waves, 
although the Z^'s contribute more to the total energy dis- 
sipation. As a matter of fact, the Z^'s develop an inertial 
range in the whole domain and the resulting cascade is more 
"aged", their dynamics is governed by the nonlinear interac- 
tions which have shorter timescales (as Z^ > Z ) and act for 
longer periods, the Z having a smaller propagation speed. 
The contribution of coherent interactions Z^Z'^^om shaping 
the energy spectra - \Zn'^ l(Akn) can be seen in the low 
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Fig. 1. — Root mean squared amplitudes of the mother wave Z+ and the reflected waves Z in km s"' (left and right panels respectively) as a function of 
heliocentric distance and time (in unit of t„.), along with the outgoing and ingoing characteristics U ±Va (red curves). Two dift'erent paths of the can be 
distinguished in the right panel: one associated to its phase speed {U - V,,, negative below the Alfvenic critical point at ai 13 Rq) and the other associated to 
reflection from the which forces the fluctuations to follow the outgoing wave (with phase speed U + Va)- 

increasing r. 

The asymptotic slope -5/3 and that one of the reflection 
dominated spectra can be deduced from the expression of 
spectral energy flux of the shell model (oi = 11/24, 02 - 
l-ai, fl3 = -15/24), 
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Fig. 2. — Energy spectra (solid and dashed lines respectively) as a func- 
tion of the perpendicular wavenumber in four different planes, as indicated 
at the bottom left of each panel (r = R/Rq). Spectra are averaged in time 
and normalized to » 1.5 X 10"^ cm^ 

power-law with the indicated slopes, plotted for reference. 
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corona, before 1.2 Rq (left panels of Figure |2]l. In fact at 
1.09 Rq, assuming waves of period T = 1000 s, e ^ 0.07 
i^max - 0.16 at the coronal base), hence in presence of lower 
frequency fluctuations (injected by the forcing) reflection is 
relatively high. At about 1 .5 Rq, Va has a maximum and re- 
flection vanishes, so that the Z are not produced locally but 
rather propagating from above. This inhibits the spectral en- 
ergy flux of that becomes negative at large scales. The 
inverse spectral transfer is responsible of the flatter part of 
the spectrum one finds at small k±. Further out reflec- 
tion is negligible and the spectra are the remnant of those pro- 
duced in the inner layers, showing a slow evolution toward 
the asymptotic state in which > 3E . The spectral energy 
flux is given by the interctions Z^Z^7^^^^ that are subject to the 
"Alfven effect" (a longer cascade timescale due to the nonlin- 
ear interactions between wavepackets propagating in opposite 
directions), producing the same spectral slopes -5/3 for E~ 
(right panels of Figure |2]i. 

Note that the largest perpendicular scale is proportional to the 
flux tube width, Ao{r) - Ao^q ^[Mr), hence the same number 
of shells spans a interval that shifts to smaller values with 



n; 



= -im [k„z^ (fliz;_^iZ;_^2 + fl2Z;_,iZ;_^2) 

+ ^„_iZ„\, [aiZUZ- + a3Z„_,z:)] 



(2) 

When reflection is negligible one can assume that the Z^ 
and Z~ are uncorrected {Z'^Z^*}t ~ 0. Assuming a power 
law for the spectral energies, E"^ cc {knlkaY*, substituting 
Z^ - 2(k„E^yi^ in eq. |2] one finds that 11^ is indepen- 
dent of the shell index n when = = -5/3. For a 
given plane, it implies a constant normalized cross helicity 
(Tc - {E^ - E^)I{E^ + E^) in the inertial range, with a cross 
helicity spectrum - E^ - E^ of the same slope as E'^, in 
contrast to EDQNM closure models, which, including non- 
local interactions in Fourier space, pred ict a distinct steeper 
spectrum, oc A;^^ dGrappin et al.lll98"3b . From the evolution 
of the spectra (right panels) one can see that the normalized 
cross helicity decreases with r because of two factors, the 
general increase of the total energy and the decrease of He at 
all scales. The former is due to the approximate conservation 
of the total wave action density. The later results from 
the competition of the linear coupling (reflection), which 
forces Z < Z^ for low frequency fluctuations, with the 
nonlinear coupling, which damps the non dominant wave 
population, the Z dDmitruk et alj 12001) . Indeed the E"^ 
spectra look the same at small wavenumbers (correspond- 
ing to n = 0, 1 , 2) in which low frequency fluctuations reside. 

When reflection is high the two fields are strongly coupled 
and one can assume that the growth of Z,7 has contributions 
both fromt the nonlinear cascade and direct generation 
through reflection x e„Z^. The reflected contribution shares 
the phase properties of Z^ so interactions due to reflection 
are coherent below 1 .2 Rq and the cascade may be isotropic 
for small wavenumbers in the kj^ - k\\ plane. For the first 7 
shells (2 orders of magnitude), one would find w„ - k„Va, 
e„ oc /I " and the following scaUng for the spectra E'^ oc ^ 
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Fig. 3. — Time averaged heating per unit mass as a function of distance for 
different root mean squared amplitudes at the coronal base Zg . Also plotted 
in red is the heating function necessary to sustain the background specified 
solar wind. H is in units of k 3 X 10"'cm^ s"'. 

oc (IVelli et al.lll99C ft. If the nonlinear cascade is some 
how inhibited, keeping fixed the total energy in the parallel 
wavenumbers, low frequencies are confined to smaller « and 
e„ oc A^"'' with p < I, yielding a steeper (flatter) spectrum 
for (£■ ) compared to the isotropic case. In the present 
simulations, ei.09Ro ^^"'^ for n < 8 which yields = 1.2, 
p^ = 2.6 as in the bottom left panel of Figure|2] 

The resulting turbulent heating is computed from 
eq. [U multiplied by (Z;f)* and summed over the shell 
index. The total energy dissipation per unit mass 
// = g/P = l/2v2;„fe,i(|Z„+|2 + |Z,7|2) increases in the 
low corona and decreases exponentially in the sub-Alfvenic 
region of the wind, as expected. For 6uq » 49 km s"' (solid 
line in Figure[3]l, it is also very close to the heating required to 
sustain the specified background wind (red line in Figure O: 
the location and height of the peak coincide, although more 
energy is dissipated in the upper corona. This would produce 
a faster wind, not altering too much the mass flux, since the 

S ieak of the dissipati on is close to the sonic critical point 
Hansteen et anil999h . 

For smaller values of Suq a better agreement is found in the 
decreasing part, the rapid increase in the low corona is still 
reproduced but the peak intensity is not attained. Despite the 
fact that the spectra and the spectral fluxes possess the same 
properties, for 6uq x 45 and 42 km s"' the peak intensity 
decreases by a factor 1/4 and 1/2 respectively (dashed and 
dotted line in Figure |3]l. It turns out that the amount of dis- 
sipation is very sensitive to the level of velocity fluctuations 
at the base of the corona. In fact, 6u is determined by the 
injection of at the coronal base but also by the response 
of the atmosphere and by the nonlinear interactions (the level 
of Z ). On the countrary the peak dissipation seems to scale 
linearly with the rms amplitude Zg although further studies 
(dependence on the frequency, on the nonlinear interactions 
in the shell model, on the imposed wind) are necessary to 
define the scaling precisely. 

4. DISCUSSION 

We have studied the propagation, reflection and nonlinear 
interaction of Alfven waves from the base of the corona up 
to 17 solar radii, well beyond the Alfvenic critical point. For 
the first time 2D shell models have been applied to account 
for nonlinear interactions in magnetically open regions on the 



sun, such as coronal holes. Thanks to such a simplification, 
compared to MHD or Reduced MHD direct numerical simu- 
lations, it is possible to follow the development of a turbulent 
spectrum in the expanding solar wind, where waves are con- 
tinuously reflected by the gradients in mean fields. Reflected 
waves are made of two components, one propagating with 
the characteristic phase speed U - Va (Z^^^^^) and the other 
following the path of the outgoing wave with spe ed U + Vg 
(Z~ ), a confirmation of prev ious Unear results (IVelli et al.l 
ll989tlHollweg & Isenberg|2007 ) which hold in a similar way 
also in the nonUnear regime. For typical coronal parameters 
we find Z - e(r, a))Z^ oc V^Z^, in contrast to Z oc AqV^ 
independet of Z'^) found in a strong turbulence regime 
Dmitruk et al.l |2002|) in which Z , Aq — > 0. Differences 
arise because the above limits impose a time scale ordering 
T'^i « T'^^i < Ts < Tcr which is not satisfied in our simulation, 
basically because Z" even in case in which the outer 
scale of turbulence does not follow the flux tube expansion 
(similar to the limit /Iq — » 0). 

Close to the coronal base one can distinguish the contribution 
to the spectral slope of the coherent nonlinear interactions 
ZanoinZ^^ which give E'^ oc fc"'^, and the one from the 
incoherent nonlinear interaction Z^^^^^Z^, giving E"^ oc k'^^^ 
(in this shell model, which includes nonlinear interactions 
only locally in Fourier space). The resulting spectra change 
with distance, starting from a coherent-interaction domi- 
nated spectrum at the coronal base and evolving toward the 
asymptotic Kolmogorov spectra at greater distance, where 
reflection is negligible. According to this model, outside the 
Alfvenic critical point, the turbulent spectra have already 
lost any feature acquired in the low corona. Note that this 
is referred to perpendicular wavenumber spectra and not to 
the frequency spectra, which on the countrary are almost 
unchanged, since their evolution is limited to the first solar 
radius above the coronal base. 

Turbulent dissipation is remarkably high in the low corona. 
Depending on the injected energy an almost complete or 
partial matching is found with the "theoretical" heating, that 
is the one required to form the imposed background solar 
wind. The best agreement is found for 6uq » 50 km s~' 
which is at the limit of observational constraints ( Chae et all 
119981) . Nonetheless, even for more conservative values 
6uq ss 40 km s"', turbulent dissipation accounts for half of 
the above theoretical heating, maintaining the same profile 
(i.e. a peak at the sonic point). This implies that the role of 
coherent interactions is fundamental in shaping the heating 
function and that turbulence and turbulent heating can not be 
neglected when studying the acceleration of the solar wind. 
The peak dissipation seems to scale linearly with the rms 
amplitude Zg although a proof of the precise scaling would 
require further studies. We finally observe that if the lower 
boundary is shifted to the base of the chromosphere, hence 
including the transition region, a stronger dissipation rate is 
expected to be found at the transition region and in the low 
corona. The Alfven speed is smaller below the transition 
region but its gradients are higher, increasing the amount of 
energy residing in the anomalous reflected component which, 
having more time to interact with its mother wave, might 
increase the spectral energy transfer 
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